A periodic microfluidic bubbling oscillator: 
insight into the stability of two-phase microflows 
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This letter describes a periodically oscillating microfoam flow. For constant input parameters, 
both the produced bubble volume and the flow rate vary over a factor two. We explicit the link 
between foam topology alternance and flow rate changes, and construct a retroaction model where 
bubbles still present downstream determine the volume of new bubbles, in agreement with experi- 
ment. This gives insight into the various parameters important to maintain monodispersity and at 
the same time shows a method to achieve controlled polydispersity. 

PACS numbers: 47.55.Dz, 47.60.+i, 83.50.Ha, 83.80.Iz 
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Formation and flow of bubbles in microfluidic systems 
attract an increasing attention i^, i^J , while poten- 

tial applications cover a wide range of fields varying from 
biology to biomedicine (ultrasound contrast agents) to 
chemistry (microreactors) . Recently these studies have 
been extended to the low liquid fraction case 0, |1] in 
which microfoams are formed. Microfoams are interest- 
ing because of the large number of interfaces that can 
potentially be used as a transport vector for amphiphilic 
molecules. Another positive aspect is the excellent con- 
trol over the volume of the dispersed phase (monodisper- 
sity) with a volume variation smaller than 3 % 0- This 
is an important factor e.g. for the synthesis of new micro- 
materials (solid microfoams). Recently a new set of op- 
erations called 'discrete microfluidics' y\ has showed the 
way to integrate microfoams in lab-on-a-chip systems. 

In this letter we will show the generation of a micro- 
foam with a periodical oscillation in its bubble volume. 
This is an unexpected effect in microfoam flows for which 
at constant input parameters both gas flow rate and bub- 
ble volume show a large periodic oscillation. We can ex- 
plain it by linking bubble volume to resistance to flow 
of a foam in a channel. This resistance depends on the 
bubble volumes downstream and on the foam topology 
therefore the volume of a new bubble is related to the one 
of its predecessors. One could state that the bubble for- 
mation process has a memory. This can be used for the 
generation of microbubble aggregates with a controlled 
and potentially tunable polydispersity. It will also play 
an important role in defining design criteria for systems 
in which monodispersity is important, both for foams 
as for emulsions with a large viscosity contrast between 
continuous and dispersed phase. The theoretical model 
provides an effective tool to predict the sensibility of the 
system to perturbations in monodispersity. It will help to 
select beforehand the right geometrical and experimental 
parameters. 

We use a conventional fiow- focusing geometry Q, |^ |^ : 
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an inlet channel for the liquid, another one for the gas, 
both ending in a small orifice followed by a straight chan- 
nel section (orifice width Wor = 100 fim, channel width 
w = I mm, height h = 100 /im, outlet channel length 16 
mm). The walls of the system are made in PDMS glued 
to a glass cover slide. A syringe pump (11 PicoPlus, 
Harvard Apparatus) is used to push the liquid (deion- 
ized water with 10% commercial diswashing detergent 
Dreft, Proctor&Gamble). For the dispersed gas phase 
we use nitrogen supplied from a pressurized tank via a 
pressure reduction valve. The exit is at atmospheric pres- 
sure. Therefore we have access to the imposed pressure 
drop over the total system (orifice and outlet channel). 
We use a CMOS camera (F131B, Alhed Vision Technolo- 
gies) connected to a macro lens to capture still images 
and to record movies of the two-phase flow. Via image 
analysis, we extract the bubble volume Vf,, bubble for- 
mation frequency /, from which we deduce the gas flow 
rate Qg^Vb f. 

Here we present a case in which for constant input pa- 
rameters Pg and Qi the bubble volume and gas flow rate 
periodically oscillate over more than a factor 2 (see Fig. 

A typical cycle starts with a foam consisting of two 
rows of bubbles (that we name F2 foam). See Fig. fflb)- 
The flow rate grows, and at the same time the bubble vol- 
ume Vf, of newly formed bubbles at the orifice increases. 
This increase is followed by a transition of the foam topol- 
ogy in the channel from two rows to one row of bubbles 
(F2 to Fl foam) (see Fig. ^c)). Now the foam starts 
slowing down while the bubble volume at the entrance 
is shrinking. This is followed by another transition, this 
time from Fl to F2 (Fig. [T]d)). When the last Fl bub- 
bles are evacuated from the channel a new cycle can start. 
Over the whole cycle the bubble formation frequency / 
stays constant. This is consistent with earlier measure- 
ments that showed that for bubble formation at high 
dispersed phase volume fraction, / only depends on the 
liquid fiow rate Qi. Therefore there is a linear relation 
between the volume of a newly formed bubble and the 
velocity v at which the foam flows: Vf, = Qg/ f — vS/f. 
Velocity is homogeneous over the channel. Image analy- 
sis shows no compressibility effects in agreement with the 
small pressures (typically a few kPa) . The foam behaves 
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FIG. 1: A periodically oscillating microfoam generator, a) 
Gas flow rate Qg (circles) (/il/s) oscillates over a factor 2.8 in 
time t (non-dimensionalised using the bubbling frequency /). 
Produced foam topology oscillates between two states: F2 (b, 
whole channel, c, near entrance) and Fl (d). The number Nfi 
of Fl bubbles (* in a)) is inversely correlated with Qg. The 



liquid flow rate Qi 
kPa. 



6.7 10 ^1/s and gas pressure Pg — 2.90 



like a plug flow. 

We stress that the bubble formation process itself is 
stable while the flow is not. We are in the high-pressure 
regime where the gas- water interface pinches off but stays 
in the orifice Non-linearities due to rapid gas- water 
interface retraction upstream of the orifice after bubble 
pinch-off like in do not play a role. Pinch-off is stable, 
both in periodicity and spatial position. For a slightly 
lower driving pressure Pg , we observe a monodisperse F2- 
foam. A slightly higher pressure leads to a monodisperse 
Fl-foam. 

A key element for the understanding of the oscillation 
are the topological transitions (from Fl to F2 and vice 
versa). The dominant force in a foam is surface tension, 
tending to minimize bubble surface for a given bubble 
volume, ft is known that a foam is usually stuck in a local 
energy minimum an does not reach the global minimum 
[Tl|. The method used for the preparation of the foam 
is the determinant factor. This also holds for our case 
in which for the same Vf, we can both have an F2 and 
an Fl foam. Furthermore the very elongated Fl bubbles 



are obviously far from minimising their surface to volume 
ratio. The formation process in this confined geometry 
plays an important role. 

We explain the topology selection mechanism by the 
interplay between bubble formation at the orifice and the 
shape assumed by the preceding bubbles. We will dis- 
cuss first Fl and F2 formation, before describing transi- 
tions. In the case of the bubble formation for an Fl foam, 
the new bubble will first form a circular shape centered 
around the orifice (see pictures in Fig. [21a). A three- fold 
vertex (three meeting films) is formed at the point where 
the wall, the new bubble and the preceeding one meet. 
These vertices slide over the channel walls. The F2 for- 
mation is a bit more complex (see Fig. |21b). Here again 
formation of a bubble starts with a spherical cap, with 
a vertex on each side. Because of the inherent asymme- 
try of an F2-foam the distance between these two ver- 
tices and the preceeding ones is unequal. Between the 
two vertices at shortest distance a Tl-transition (bubble 
neighbour switching |l2l|) takes place alowing the newly 
formed foam to relax to the F2 state. 

The transition from Fl to F2 takes place when the 
volume Vb of the new bubble becomes that large that the 
inter-vertex distance is too large to allow a Tl (see Fig. 12 
c). We will call this volume V12: the threshold volume for 
transition from Fl to F2. The opposite transition takes 
place at a much smaller bubble volume because of the 
absence of preexistent asymmetry: bubble volume has 
to reduce a lot (smaller than V21 ) before the inter- vertex 
distance is that small that the situation becomes unstable 
(probably because of interaction with previous film), that 
a top or bottom vertex slides faster, triggers a Tl upon 
contact with the previous film, thereby completing the 
transition ^d). 

For our channel geometry the transition Fl to F2 takes 
place at Vi2/hw'^ — 0.193 and the inverse at V2i/hw'^ — 
0.536. There is a large hysteresis (a factor 2.8 between 
these two volumes). The difference in transition bubble 
volume probably depends heavily on the liquid fraction 
as this is the determinant factor for the vertex distance at 
which a Tl rearrangement is triggered Therefore 
we expect that the amplitude of the oscillation will be less 
pronounced for wetter foams. 

We complete the explanation of the oscillating behav- 
ior by making a link between foam topology and rheology. 
Driving pressure Pg is related to foam velocity v by the 
channel resistance to flow (dissipation). The foam speed 
is oscillating for constant driving pressure, therefore the 
channel resistance must oscillate as well. The principal 
source of dissipation for foam flow in a channel is the 
sliding of the liquid films between bubbles over the chan- 
nel walls [14L Il5| . Channel resistance scales linearly with 
the total length of these films in the channel. Orienta- 
tion plays an important role. For a film of length L the 
relevant length scale is the projected length = L cos a 
with a the angle between the normal vector of the film 
and the flow direction. Introducing the capillary number 
Ca = fxv/a containing the bubble velocity v (estimated 
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a) Fl bubble formation 





b) F2 bubble formation 




c) Transition F2 to Fl 




d) Transition Fl to F2 



with the following dimensionless quantities: 14 = H/ hw^ 
the bubble volume, Rtotai = J2k=n-Nt L^^^'' /w the to- 
tal resistance of the films, Ch/ = iiwf/cr and P = 
^Pchanneih / 2a the experimental parameters related to 
frequency and pressure. The sum of the projected lengths 
Rtotai, takes into account the Nb bubbles in the channel. 
In addition, we enforce foam topology transitions at the 
bubble volumes V12 and V21 , which completes the model. 

Essentially we calculate Vi,^"^^'' by numerical integration 
over all the bubbles in the channel. Therefore a bubble 
influences the volume of its successors during its presence 
in the channel, so while travelling from the orifice to the 
channel exit. Thereby creating a retroactive effect. 

A comparison between experiment and model is shown 
in Fig. 13 confirming that we capture the essential of the 
physics. The model reproduces the oscillation frequency, 
the gas flow rate and number of bubbles in the channel 
nicely. 



FIG. 2: Formation of the foam topology Fl (a) and F2 (b), 
and tiie transitions F2 to Fl (c) and Fl to F2 (d). 



as w ~ Qg/S) and liquid viscosity /i, the pressure drop 



writes 



^Pchar, 



(1) 



with S = hw the channel cross section area, and A a 
numerical constant 0, 0, 0| . 

Using geometrical considerations we can calculate the 
film length per bubble, projected on the direction normal 
to flow, for the two distinct topologies: 
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taking into account the bottom and top wall of the chan- 
nel, and neglecting the side walls. The largest possi- 
ble bubble volume in the F2 topology is Vb/hw^ — \/3. 
Larger volumes would violate Plateau's law saying that 
three films must always meet at equal angles. Therefore 
for all possible bubble areas L^pi > L^p2- 

Summing over the projected lengths of all the bubbles 
in the channel gives the total resistance prefactor and 
thereby the capillary number Ca (and gas flow rate Qg) 
at a given applied pressure drop, from equation 

To model the foam flow we proceed in the following 
way. The bubbling frequency / being constant whatever 
the bubble size, we predict the volume of a new bubble 
to be Vb = Qg/ f- Using equation we obtain the gas 
flow rate Qg as a function of the bubbles downstream in 
the channel, and thus the next bubble volume, at a given 
discretised time tf — n, from 
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FIG. 3: a) Comparison between experimental data (circles) 
and the model (continuous line) for the bubble volume VJ, 
at the entrance of the channel. Time is normalised by the 
bubbling frequency /. b) Bubble volume against the total 
number of bubbles in the channel A^;,. (o): model, (o): ex- 
periment. Filled symbols indicate that newly formed bubbles 
form an Fl topology, open symbols the generation of F2. Pa- 
rameters A = 67, CaJ = 8.9 10"", P = 8.2, Vi2/hw^ = 0.25, 
V21 /hw'^ = 0.45. 
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